Experimental and Atomistic Theoretical Study of Degree of Polarization 
from Multi-layer InAs/GaAs Quantum Dot Stacks 



Muhammad Usman, 1 '^] Tomoya Inoue, 2 Yukihiro Harda, 2 Gerhard Klimeck, 3 and Takashi Kita 2 

1 Tyndall National Institute, Dyke Parade, Cork, Ireland 
2 Department of Electrical and Electronic Engineering, 
Graduate School of Engineering, 1-1 Rokkodai, Nada, Kobe 657-8501, Japan 
' Purdue University, Network for Computational Nanotechnology, 
Birck Nanotechnology Center, 1205 W. State Street, West Lafayette, IN 47907, USA 

Recent experimental measurements, without any theoretical guidance, showed that isotropic po- 
larization response can be achieved by increasing the number of QD layers in a QD stack. Here 
we analyse the polarization response of multi-layer quantum dot stacks containing up to nine quan- 
tum dot layers by linearly polarized PL measurements and by carrying out a systematic set of 
multi-million atom simulations. The atomistic modeling and simulations allow us to include cor- 
rect symmetry properties in the calculations of the optical spectra: a factor critical to explain the 
experimental evidence. The values of the degree of polarization (DOP) calculated from our model 
follows the trends of the experimental data. We also present detailed physical insight by examining 
strain profiles, band edges diagrams and wave function plots. Multi-directional PL measurements 
and calculations of the DOP reveal a unique property of InAs quantum dot stacks that the TE 
response is anisotropic in the plane of the stacks. Therefore a single value of the DOP is not suf- 
ficient to fully characterize the polarization response. We explain this anisotropy of the TE-modes 
by orientation of hole wave functions along the [110] direction. Our results provide a new insight 
that isotropic polarization response measured in the experimental PL spectra is due to two factors: 
(i) TMooi-mode contributions increase due to enhanced intermixing of HH and LH bands, and (ii) 
TEno-mode contributions reduce significantly due to hole wave function alignment along the [IlO] 
direction. We also present optical spectra for various geometry configurations of quantum dot stacks 
to provide a guide to experimentalists for the design of multi-layer QD stacks for optical devices. 
Our results predict that the QD stacks with identical layers will exhibit lower values of the DOP 
than the stacks with non-identical layers. 

PACS numbers: 78.67.Hc, 73.22.Dj 



Introduction and Problem Background: The de- 
sign of optical devices such as semiconductor optical am- 
plifiers requires polarization insensitive optical emissions. 
The polarization response of quantum dots is measured 
in terms of degree of polarization (DOP), defined as: 
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DOP = JLJi (1) 

TEj_- growt h + TM\\_ growth 

Here TEj__ srotot ^ refers to traverse electric mode in a 
direction perpendicular to the growth direction ([001]) 
and TM||_ growt h refers to traverse magnetic mode along 
the growth direction ([001]). The value of DOP depends 
on the chosen direction for the TE-mode in the plane 
of quantum dot. While the previous studies of polariza- 
tion response of quantum dots provide only one value of 
DOP corresp onding to a chosen direction for the TE- 
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we associate DOP with the direction of 
TE-mode, for example TE 110 — > DOPn and TE Tl0 — > 
DOPj 10 , and show that the DOP depends highly on the 
chosen direction. Only one DOP value is not sufficient to 
fully characterize the polarization response of quantum 
dots systems. 

The InAs quantum dots obtained from the Stranski- 
Krastanov self-assembly growth process typically have 
a flat shape i.e. base diameter is typically 4-5 times 



larger than the height. In such quantum dots, com- 
pressive biaxial strain splits the heavy hole HH and 
light hole LH bands by more than lOOmeV. As a re- 
sult, only the TE-mode can couple and the TM-mode 
is very wealf^. The polarization response of such sys- 
tems is highly anisotropic, TE-mode 3> TM-mode and 
DOP— >1.0. To achieve the desired isotropic response 
(DOP=0) for the design of optical devices, significant 
tuning of QD geometry, band structure manipulation, 
and/or strain engineering are required. 

During the last few years, several techniques have been 
explored to achieve polarization insensitive optical emis- 
sion from InAs quantum dot (QD) samples. These meth- 
ods include overgrowing the InAs QD samp les by an In- 
GaAs strain relaxing capping layer (SRCLjEE growing 
large stacks of QDs in the form of colu mnar Q Da^"^, bi- 
layer^, tri-layer QEP^, multi-layer stacks^^, and band 
gap engineering by including dilute nitrogen ]\PS, phos- 
phorous P^3, and antimony SfP3 impurities. T. Kita et 
aZP and T. Saito et demonstrated that an isotropic 
polarization response can be obtained by growing colum- 
nar QDs consisting of nine QD layers. 

Recent experiments by T. Inoue et a/! 12 * 19 * 36 ^ showed, 
without any theoretical guidance, that similar tuning of 
polarization properties is possible in regular InAs QD 
stacks where the QD layers are geometrically separated 
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FIG. 1: (a)TEM images of three, six, and nine QD layer stacks grown by solid-source molecular-beam epitaxy^. The layers of 
InAs QDs separated by GaAs intermediate layers are stacked with clear wetting layer interfaces, (b, c, d) Linearly polarized 
PL measurements for TEno, TEt 10 , and TMooi-modes for single QD layer, (e, f, g) Linearly polarized PL measurements for 
TEno-, TEjio-, and TMooi-modes for nine QD layers sample. The PL plots in the figures (b) to (g) are from independent 
experimental measurements. Each one of them is normalized to its highest peak, (h) The schematic diagram of the simulated 
system. InAs quantum dot layers are embedded inside GaAs buffer. Each layer consists of a dome-shaped quantum dot on top 
of 1.0ML InAs wetting layer. QD layers are denoted by Ljv, where N being the number of QD layers in the stack. The wetting 
layers are separated by 4.5nm GaAs buffer. 



by thin GaAs spacers. Such multi-layer QD stacks have a 
twofold advantage over columnar QDs: (i) a moderately 
thick GaAs spacer between the QD layers allows a pre- 
cise control of overall QD shape and size, (ii) a reduced 
strain accumulation results in isotropic polarization re- 
sponse with fewer number of QD layers in the stack. The 
experimental results are provided for QD stacks contain- 
ing three, six, and nine QD layers. The results indicate 
that the DOPno takes up the values of +0.46 and -0.60 
for the samples containing six and nine QD layers, respec- 
tively. A change of sign for the DOPno implies that an 
isotropic polarization response (DOP^O) can be achieved 
by engineering the number of QD layers in the stack. 

The previous experimental measurements^ indicated 
a prospect to achieve polarization insensitive response, 
however no theoretical study is available to-date to pro- 
vide physical insight for the design of these complex 
multi-million atom nano-structures. Furthermore the ex- 
periments only analyse TEno-mode whereas PL mea- 
surements along multiple in-plane directions are required 
to fully characterize the polarization properties of these 
QD systems. 

This article aims to thoroughly investigate the optical 
spectra of multi-layer QD stacks by in-plane polarization 
measurements and atomistic theoretical analysis. Our 
PL measurements reveal unique properties of these QD 
stacks containing strongly coupled electronic states and 
the atomistic theory explains the experimental evidence 
very well. The modeling and simulations also provide 
data for DOP for different geometry configurations and 
in-plane directions for the TE-mode to explore the design 
space and provide a guide for future experiments. 



Experimental Procedure: The QD samples are 
grown on an undoped [001] GaAs subtrate using solid- 
source molecular-beam epitaxy. First a thick layer of 
GaAs buffer layer is grown at 550°C, followed by rep- 
etitions of the InAs QD layers and GaAs spacer layers 
with growth interruptions of 10s after each GaAs spacer 
layer reached a thickness of 16MLs. The nominal thick- 
ness of the InAs is kept 1.9ML. Finally the QD layers 
are capped with a lOOnm thick GaAs layer. Further 
details of our growth process can be found in earlier 
publication s] 12 * 19 ! 23 * 36 !. The crystallographic properties 
of the stacked QDs are examined using a cross-sectional 
transmission electron microscope (TEM). The TEM im- 
ages of the samples containing three, six, and nine QD 
layers are shown in Fig. [TJa^P^l 

Next, we perform linearly polarized PL measurements 
at room temperature to investigate the polarization de- 
pendent optical spectra. The laser diode excitation is at 
659nm wavelength. The detailed set up and measure- 
ment procedure is described by T. Inoue et aP^. The 
polarization dependent PL spectra are shown in Fig. [I] 
for a single QD layer (b, c, d) and nine QD layers (e, f, 
g). 

Fig. [ijb, c, d) indicate for a single quantum dot that 
TE 110 ~ and TEj 10 show a similar magnitude, and that 
TMooi has a much smaller response than both in-plane 
TE-modes. This is typical for InAs QDs and is due to the 
compressive biaxial strain that splits the HH-LH band 
edges resulting in HH-type valence band states close to 
band gap. Previous theoretical and experimental studies 
on single QDs have shown similar properties. 

The PL spectra in the Fig. [jje, f, g) on contrary reveal 
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an interesting and unique property of the InAs/GaAs 
QD stacks that the TEno response is significantly less 
than the TEj 10 response. No previous evidence exists 
for such TE-mode anisotropy in InAs/GaAs QDs, though 
Podemski et alJ^ have reported similar results for InP- 
based columnar quantum dash structures. The measured 
difference between the two TE-modes is such that even 
the DOP for the same QD stack could have different 
signs when measured along the [110] and [110] direc- 
tions. The reason for such TE-mode anisotropy is highly 
non-intuitive and requires modeling and simulations of 
these large QD structures. The modeling must include 
the correct symmetry properties because the long range 
strain effects the room temperature material properties 
and atomistic structure resolution. The multi-million 
atom simulations presented in this article explains the 
experimental measurements in terms of hole wave func- 
tion alignments along the [110] direction. 

Theoretical Model: The theoretical modeling of QD 
stacks containing up to nine QD layers posses a twofold 
challenge: first it requires an atomistic model that can 
calculate electronic and optical properties including cor- 
rect symmetry and interfaces. Secondly the large size 
of QD stacks requires calculations to be done over mil- 
lions of atoms to properly incl ude t he long-range effects 
of strain. We use NEMO 3-EPEl to analyse the elec- 
tronic and optical properties of multi-layer QD stacks. 
NEMO 3-D is based on fully atomistic calculations of 
strain and electronic structure: the strain is calculated 
by using atomistic valence force field (VFF) modef^, in- 
cluding anharmonic corrections^ to the Keating poten- 
tial and the electronic structure is calculated by solving 
a Hamiltonian inside twenty band sp 3 d 5 s* basis^. 

NEMO 3-D has been designed and optimized to be 
scalable from single CPU to large number of proces- 
sors on most advanced supercomputing clusters. Ex- 
cellent MPI based scaling to 8192 cores/CPU has been 
demonstrated^. The atomistic modeling techniques and 
parallel coding scheme implemented in NEMO 3-D al- 
lows it to simulate large QD stacks with realistic geome- 
try extension and symmetry properties. Past NEMO 3-D 
based studies of nano-structures include (i) s ingle InAs 
QD with InGaAs SRClP, InAs bilayer QDs^U, and InAs 
multi-layer QD stacks^, (ii) valley splitting in miscut Si 
quantum wells on SiGe substrate^, (iii) Stark effect of 
single P impurities in SP^l, and (iv) dilute Bi impurities 
in GaAs and GaP materialman. 

Simulated System: The theoretical analysis of the 
multi-layer QD stacks containing up to nine layers is car- 
ried out through a set of systematic simulations. The 
geometry parameters of the QD samples are extracted 
from the TEM images shown in the Fig. [lja) which in- 
dicates that all the QDs in the stacks arc of nearly iden- 
tical size. Some top most layers in the nine quantum dot 
stack appears to be relatively small and are only par- 
tially growrP^l The schematic diagram of the system 
modelled in our simulations is shown in Fig. [ljh). Mul- 
tiple QD layers separated by 4.5nm thick GaAs buffer 



are embedded inside a large GaAs matrix. This wetting 
layer-to-layer separation is experimentally optimized to 
obtain uniform vertical QD stacks from the self-assembly 
growth process^. Each QD layer consists of a dome- 
shaped InAs QD with circular base, lying on top of a 
1.0ML InAs wetting layer. The QD systems with single, 
three, six, and nine QD layers will be labelled as L 1: L 3 , 
Li6, and Lg, respectively. 

The size of QDs extracted from the TEM images in- 
dicates a base diameter of ^20nm and height of ~4nm. 
Since the height of the quantum dots in the TEM images 
is not very clear, we choose to simulate three different QD 
geometries in our theoretical study: (i) All QD layers are 
identical with 20nm base diameter and 4nm height of the 
QD in each layer. All the results presented are for this 
system unless other dimensions are specified, (ii) All QD 
layers are identical with 20nm base diameter and 3.5nm 
height of the QD in each layer and (iii) the size of QDs 
increases from the lower to the upper layers: the base 
diameter increases by lnm and the height increases by 
0.25nm. We choose this last syste m (hi ) because past 
experiments of multi-layer QD stacks^GI] have indicated 
an increasing size of QDs when multi-layer QD stacks are 
grown by the self-assembly process, so it is interesting to 
theoretical investigate this geometry. We also simulate 
a system L g in which the height of QDs in each layer 
is 4.5nm. This would mean that the top of each quan- 
tum dot will be touching its upper adjacent wetting layer, 
approaching columnar QD limit. The theoretical results 
presented here for various QD geometries provide a guide 
for the experimentalists to understand the dependence 
of DOP on the QD geometry since experimental inves- 
tigations of QD stacks for polarization response are still 
under way. 

The QD layers are embedded inside a sufficiently large 
GaAs buffer to ensure proper relaxation of atoms and 
to accommodate the long-range effects of strain. The 
size of the largest GaAs buffer for the system containing 
nine QD layers Lg is 60x60xl06nm 3 , consisting of ~24.4 
million atoms. Mixed boundary conditions are applied 
in the strain minimization: the substrate is fixed at the 
bottom, the GaAs matrix is periodic in lateral dimen- 
sions, and the capping layer is free to relax from the 
top. The electronic structure calculations use a separate 
subdomain with closed boundary conditions to reduce 
the computational burderP^, with a surface passivation^ 
which avoids artificial surface states in the atomistic rep- 
resentation. 

Electron Wave functions form molecular 
states: Figure [2] shows the plots of the lowest conduc- 
tion band state Ei for all of the four QD systems under 
study. From the top views of the wave functions (sec- 
ond row), it is evident that the lowest electron state is 
of s-type symmetry. The side views of the wave func- 
tions (first row) show that the electron state forms a 
hybridized (molecular) state in L3, Lg, and Lg stacks, 
and is spread over all of the quantum dots. This is due 
to strong coupling of quantum dots at 4.5nm separation. 
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FIG. 2: Plots of the lowest conduction band state, Ei, for 
quantum dot systems Li, L3, Lg, and Lg. First row: the side 
view of the plots is shown. Second Row: the top view of the 
plots is shown. The intensity of the color in the plots indicates 
the magnitude of the wave function: the red color represents 
the highest magnitude and the light blue color represents the 
lowest magnitude. The dotted circles are marked to guide the 
eye and indicate the boundary of the base of each QD. 



The presence of the s-like electron wave function in all 
of the quantum dots implies that only the details of hole 
wave functions inside the quantum dot stack will deter- 
mine the optical activity of a particular quantum dot in- 
side the stack. This is different from the previous study 
of bilayers^l where a weak coupling of quantum dots at 
lOnm separation resulted in atomic like electron wave 
functions. In the bilayer system, both the electron and 
the hole wave functions determined the optical activity 
of the quantum dots. Previous studies on the identical 
bilayepSD and stacks with seven identical QDsP^ showed 
that the strain tends to push the electron states towards 
the lower QDs in such systems. The electron wave func- 
tions for L3 and Lg systems follow this trend. Here we 
find that this trend is no longer true for stacks with nine 
QD layers, Lg, where the electron wave function Ei van- 
ishes around the edges of the stack due to the larger strain 
magnitude there (see Fig. [2] first row for Lg). 

Hole wave functions exhibit atomic charac- 
ter: Although the electrons (lighter-mass particles) 
are strongly influenced by the inter-dot electronic and 
strain couplings of QDs and exhibit tunnelling across the 
quantum dots forming molecular like hybridized states, 
the holes due to their heavier mass remain well con- 
fined inside the individual dots and do not show any 
hybridization^. For example, Fig. |3]show the side views 
of the lowest five valence band wave functions in the Lg 
system. The horizontal dotted lines are plotted to mark 
the positions of the base of the QDs and helps to de- 
termine the location of a particular hole wave function 
inside the stack. In this system, Hi and H3 are inside 
QD 2 , H 2 and H 5 are in QD 3 , and H 4 is in QD 8 . The 
location of a hole state inside a QD stack is relatively 



hard to determine and is strongly influenced by geome- 
try of the QD stack i.e. QD base diameter, QD height, 
QD layer separation etc. Ultimately the strain profile 
that controls the strength of the coupling between the 
QD layers inside the stack determines the position of the 
hole states inside the stacks. 

Hydrostatic and Biaxial Strains: Figure [4] plots 
the hydrostatic €h=€xx + ^yy + <= Z z (dotted lines) and 
biaxial strain Gb=^xx + &yy —2 e zz (solid lines) profiles 
along the [001] direction through the center of the quan- 
tum dots. The hydrostatic strain exhibit a very slight 
change from Li to Lg. The biaxial strain however sig- 
nificantly changes as the QD stack height increases. For 
a single QD, Li, the biaxial strain is highly negative in- 
side the QD region. As the vertical size of the QD stack 
is increased by adding QD layers, the biaxial strain at 
the center of the stack reduces. For the Lg stack in the 
Fig. Qd), the biaxial strain at the center of the stack 
approaches zero. The reason for such behaviour of the 
biaxial strain is that in general the InAs unit cells inside 
the QD region tend to fit over the GaAs matrix by an 
in-plane compression and an elongation along the [001] 
direction. This results in highly negative biaxial strain 
as can be observed for a sing le QD in the Fig.ga). How- 
ever, when the size of the stack increases, the unit cells 
of InAs around the center of the stack feel lesser and 
lesser compressive force from the surrounding GaAs. As 
a result, the vertical lattice constant of the InAs starts 
matching with the GaAs and hence the biaxial strain 
tends to change its sign around the middle of the QD 
stack. Similar strain profiles were calculated in an earlier 
study about columnar QDs by T. Saito et alW. 

Increased HH-LH Mixing: The minor change in 
the magnitude of the hydrostatic strain (as Li — > Lg) 
implies that the lowest conduction band edge will expe- 
rience very small change as they are only effected by the 
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FIG. 3: Plots of the highest five valence band state, Hi, H2, 
H3, H4 and H5, for quantum dot systems L9. Only the side 
view of the plots is shown. The horizontal dotted lines are 
marked to guide the eye and indicate the base of the QD layers 
in the stack. The intensity of the color in the plots indicates 
the magnitude of the wave function: the red color represents 
the highest magnitude and the light blue color represents the 
lowest magnitude. The energies of the valence band states and 
the differences between the energies of the adjacent levels are 
also mentioned. 
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FIG. 4: The plots of hydrostatic (6h= 
biaxial strain (€b=€xx + € yv —2 G z 
the center of the quantum dot system along the [001] direc- 
tion. The hydrostatic component is dominantly negative in- 
side the QD indicating strong compression of the InAs and 
alomost zero outside the QD. As the QD stack height in- 
creases, the biaxial strain in the QD evolves from negative 
to zero with a small increase in positive contributions in the 
capping layer. 



hydrostatic component^. The valence band edges are af- 
fected by both the hydrostatic strain as well as the biaxial 
strain. The impact of strain on the highest two valence 
band edges, HH and LH, is analytically expressed as: 



SE LH = a v £ H 



2 



(2) 



(3) 



Here a v and b v are deformation potential constants 
for HH and LH band edges. The values for these con- 



-0.4 | 

10 15 20 25 30 35 40 45 
0.6r 



50 55 



> 

>. 

aj 
c 



0.4 
0.2 


■0.2 
-0.4 



(b)Ls 
LH 



.... 



HH 



HH 



LH 



T* 



LH 



HH 



3 



> 

<_ 
C 



10 15 20 25 30 35 40 45 50 55 
0.6 
0.4 
0.2 


-0.2 
-0.4 
1 



:(c)L 6 








•■~L. ' i 
iii .•* • 


HH 

»»>«*»■ •« 

LH 


yS: . 


:£"\lh 


- Hir-^ \ * 






- % \ 


% 


% \ 


X HH 



10 



20 



30 



40 



50 



60 



70 



; 

C 



0.6 
0.4 
0.2 


0.2 



-0.4 



(d)L 9 



i^^i \ }f — * * 

, US SC XC S*C XZ 5«C 




10 20 30 40 50 60 70 80 90 

Distance along [001] (nm) 

FIG. 5: Plots of local band edges for the highest two valence 
bands (HH and LH) through the center of the QD along the 
[001] direction for Li, L3, Lg, and L9. The reduction in the 
magnitude of biaxial strain results in larger HH/LH intermix- 
ing (lesser separation) as the size of the stack increases. 



stants for InAs systems are a v = l.OeV and b v = -1.8eV, 
respectively 1 . From the equations 2 and 3, it is evident 
that the magnitude of the gg determines the HH-LH 
splitting. For a single QD layer, due to large negative 
value of £b, the HH and LH band edges will be con- 
siderably separated inside the QD region. This will in- 
duce dominant HH character in the highest few valence 
band states which will be closer to the HH band edge. 
As the magnitude of Gb decreases, the HH-LH splitting 
reduces, increasing LH component in the valence band 
states. For the Lg system, the nearly zero magnitude of 
the G b around the center of the stack implies that the 
HH and LH bands will be nearly degenerate around the 
center of the QD stack. The valence band states will 
therefore be of highly mixed character, consisting of con- 
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tributions from both the HH and the LH bands. 

Figure |5ja-d) plots the highest two local valence band 
edges, HH and LH, for all of the QD systems under study 
along the [001] direction through the center of the QDs. 
Highly negative biaxial strain in Li results in ~152meV 
splitting of HH and LH bands within the QD region. 
As the biaxial strain around the center of QD stack de- 
creases (approaching towards zero), the HH-LH splitting 
around the center of the stack also decreases to ~74meV, 
~32meV, and <28meV for the L 3 , L 6 , and Lg systems, 
respectively. 

The HH and LH character of a particular valence band 
state in the tight binding formulation can be estimated 
as follows: If the amplitudes of the p^, p y , and p z orbitals 
at any atomic site are & x<u /d, & y ,u/di and a z ,«/d respec- 
tively (where the subscripts u and d refer to up and down 
spin) , then the HH contribution is approximately propor- 
tional to \a XtU —ia VjU \ +\a Xt d+ia Vt d\ 2 summed over all the 
atomic sites. The LH contribution is approximately pro- 
portional to |a 2jU | 2 -(-|a Zi d| 2 summed over all the atomic 
sites. By using these expressions, we estimate that the 
jjY ratio for the highest valence band state (Hi) de- 
creases from ^108 for the Li to ~15.8, ~12.3, and ~10.6 
for the L3, Lg, and Lg systems, respectively. This clearly 
points towards an increasing LH character in the valence 
band states as the height of the QD stack increases. 

HH/LH intermixing implies TM 00 i-mode in- 
creases: In a QD system, the HH states consist of con- 
tributions from p x and j> y orbitals and the LH states 
consist of contributions from p x , p y , and p z orbitals^. 
These configurations imply that the TM-mode (which is 
along the z-direction) will only couple to the LH states. 
The large splitting of the HH-LH bands (see Fig. [5]ja)) 
resulting in a weak LH contribution in the Li system will 
result in very weak TMooi-mode for this system. Thus 
from Eq. (1), the DOP will be nearly 1.0 and the polar- 
ization response will be highly anisotropic. As the size of 
the QD stack increases, the larger intermixing of HH and 
LH bands increases LH contribution in the valence band 
states. This will result in increase of TMooi-mode of op- 
tical transitions reducing the anisotropy of DOP bringing 
it closer to 0. 

Optical intensity functions, f(\): Figure [6] plots 
the optical intensity functions computed from our model 
as a function of the optical wavelength for various quan- 
tum dot systems. The calculation of the optical intensity 
function is done as follows: first we calculate optical tran- 
sition strengths by using Fermi's golden rule^for TEno-, 
TEj 10 -, and TM oi-modes between the lowest conduction 
band state Ei and the highest five valence band states Hi , 
H2, H3, H4, and H5 for a particular QD system. Next 
each optical transition strength is artificially broadened 
by multiplication with a Gaussi an dis tribution centred at 
the wavelength of the transitio n 1 38 ! 39 !. Finally we add all 
of the five Gaussian functions to calculate the total opti- 
cal intensity function, f(X ). The complete expression for 
the optical intensity function, plotted in the Fig. [6]ja-k), 
is given by equations 4 and 5: 



/(AWH,=£(T^)-e~ C ^)'" (4) 

i=l 

where, 

T El = (TEft - H < /TEg- H * /TM*\- H < ) (5) 



The examination of the optical intensity plots in the 
Fig. [6] reveals that the TMooi-mode indeed increases 
as the size of QD stack is increased: Li — > Lg. This 
is in general true for all of the geometry configurations 
considered and is a direct consequence of the change in 
the biaxial strain component that increases HH and LH 
intermixing as discussed earlier. 

Increase in TMqqi only partially contributes to- 
wards isotropic polarization: Fig. [6] shows that the 
increase in the TMooi-mode only partially helps towards 
an isotropic polarization response. This is in contrast 
to a general notion where it is described that the in- 
crease in the TMooi-mode is mainly responsible for the 
isotropic polarization. The reason for such understand- 
ing is that th e previous theoreticaP or the experimental 
studies^ 10 ! 12 ! of the DOP have assumed only one direc- 
tion for the TE-mode. However, our PL measurements 
in the Fig. flje-g) show that the TE-modes along the 
[110] and [110] have significant anisotropy in the plane 
of the QD stack. The theoretical model shows that in 
fact a major contribution to achieve isotropic polariza- 
tion response in these systems stems from a suppressed 
TEno-mode rather than an increased TM oi-mode. Fig. 
[6] shows that irrespective of QD geometry, the increase in 
the TMooi-mode is insufficient to reverse the sign of the 

DOPno (+ -)■ 

The TE-mode is highly anisotropic in the plane of the 
QD with the magnitudes of the TEno- and TEii -modes 
becoming very different as the QD stack size increases. 
For a single QD system, Li, TEno ~ TE no and TM oi 
is very weak. Hence the measured and calculated DOP 
is highly anisotropic (close to 1.0) irrespective of the di- 
rection for the TE-mode. As the QD stack size increases, 
the TMooi-mode also increases partially contributing to 
reduction in the DOP. However, at the same time, the 
TEno-mode decreases considerably such that for Lq and 
Lg, it becomes smaller than the TMooi-mode. The reason 
for such a drastic decrease in the TEno-mode is the ori- 
entation of hole wave functions along the [110] direction 
for the Lg and Lg systems. 

Hole wave functions are oriented along [110] 
and [110] resulting in TM QO i > TEno: Figure [Tlplots 
the top views of the highest five valence band states Hi , 
H 2 , H 3 , H 4 , and H 5 for the QD systems Li, L 3 , L 6 , and 
Lg. The five hole wave functions for the Li system have 
an almost uniform distribution inside the QD region with 
nearly symmetric shape. Such symmetry will result in 
approximately equal magnitude of TE-mode along the 
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FIG. 6: The plots of optical intensity functions, 
/(A), are shown for various QD systems. The op- 
tical intensity functions in each case are computed 
from equations 4 and 5. The rows of the figure 
represent the QD system: first row=Li, second 
row=Z/3, third row=L6, and fourth row=Lg. The 
columns of the figure represent the geometry of 
the particular system. First and second columns 
show results for identical QD stacks with 3.5nm 
and 4.0nm QD heights, receptively, and 20nm 
base diameter. The third column is a special case 
where we consider increasing size of QDs for L3 
and L6 systems. For L9, we again simulate identi- 
cal QDs, but each with height 4.5nm. This is the 
case when regular QD stack approaches columnar 
QD shape. 



[110] and [110] directions, as evident in the first row of 
the Fig. [7j For the stacks with three, six, and nine QD 
layers, the distribution of the hole wave functions is ori- 
ented along the [110] or [110] directions. This [110]/[110]- 
symmetry is mainly due to the strain and piezoelectric 
potentials that lowers the overall sy mmet ry of the QD 
system and favours these two direction d 1 * 34 !. To verify the 
impact of the strain and piezoelectricity, if we conduct a 
numerical experiment and switch off their contributions 
in the electronic structure calculations, the TEj 10 /TEn 
ratio for the Lg system decreases from 10.92 to 0.85 and 
it decreases from 10.73 to 3.2 for the Lg system. Similar 
distributions of the hole wave functions are observed for 
bilayeriPand a single QD layer with aspect raticE^l (H/B) 
> 0.25. 

The orientation of the hole wave functions determines 
the magnitude of the TEnrj- and TEj 10 -modes since the 
lowest electron wave function (see Fig. [2]) is symmet- 
rically distributed. For L 3 , the hole wave functions Hi, 
H2, and H5 are oriented along the [110] direction, while 
the other two hole wave functions H3 and H4 are oriented 
along the [110] direction. The orthogonal distributions of 
the hole wave functions in L3 will result in similar [110] 



and [110] TE-modes. The cumulative summations for 
the TEno and TEi 10 -modes arising from these five hole 
states are indeed nearly equal, with TEji magnitude 
being slightly larger as can be seen in the Fig. |6jd). 

All of the highest five hole wave functions are oriented 
along the [110] direction in case of the Lq and Lg systems. 
This results in a strong reduction of the TEno-mode. 
The significant reduction in the TEno-mode turns out 
to be even smaller than the magnitude of TMooi-mode 
as can be seen in the Fig.^f, g, j). This change in relative 
magnitude results in a flip of sign (+->-) for DOPn as 
indeed measured in the experiment^. Hence we conclude 
that the isotropic polarization response demonstrated by 
the experiment is a result of two factors: (i) increase in 
the TMooi-mode due to enhanced HH-LH intermixing 
and (ii) the reduction of TEno-mode due to orientation 
of holes along the [110] direction. 

Here we want to point out that the relative magnitude 
of the TEjio-mode does not reduce as the size of the 
QD stack is increased. Even for the largest system un- 
der study, Lg, the TEj 10 -mode is much stronger than the 
TMooi-mode. That means if the experimental measure- 
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FIG. 7: Top view of the plots of the highest 
five valence band states Hi, H2, H3, H4, and 
H5 are shown for Li, L3, Lg, and Lg systems. 
The intensity of the color indicates the magni- 
tude of the hole wave functions: the red color 
indicating the largest magnitude and the light 
blue color indicating the smallest magnitude. 
The dotted circles are marked to guide the eye 
and indicate the boundary of the QD bases. 



merits are performed for DOPi 10 , they should still show 
anisotropy. This is verified by our PL measurements on 
Lg system shown in Fig. [TJf ) . These PL measurements 
on L9 system indicate a positive value for the DOPi 10 
(TEj 10 > TM001) and a negative value for the DOPno 
(TE 110 < TMqoi). 

Comparison with experimental PL data: Table [I] 
summarizes the calculated values of DOP from our model 
along different directions and compares it with the exper- 
imental PL measurements. Here we provide theoretically 
calculated values of DOP along [100] and [010] directions 
for comparison purpose because some recent experimen- 
tal studies have chosen these directions for investigation 
of the DOP and characte rizatio n of the polarization re- 
sponse of the QD systemsJffiSmi. 

The hole state symmetries shown in the Fig. [7] show 
strong alignment in the [110] and [110] directions indicat- 
ing that symmetries to the [100] and [010] directions will 
be almost equivalent. Therefore the values of the DOP 
are nearly equal for the [100] and [010] directions as men- 
tioned in the table [I] The PL measurements along these 
two directions will not exhibit isotropic polarization even 
for the Lg system. Similar results were found in an earlier 
experimental study by P. Ridha et al. 33 . In their study 
of the polarization properties of multi-layer stacks based 
on TEgio-niode, they conclude that such systems can not 
provide isotropic polarization and columnar QDs are the 
only choice. 

For the system Li containing only a single QD layer, 
our calculated DOP values in all the directions are 
very close to 1.0, exhibiting a very strong polarization 
anisotropy (TE-mode 3> TM-mode). The experimen- 
tally measured values for this system are relatively low 
(^0.7), but they also show that the value of the DOP 
remains nearly the same irrespective of the measurement 



direction. While the reason for difference between our 
calculated value and the measured value for this QD sys- 
tem is not v ery clear, other theoretical studies using k.p 
method^ on single QD layers with similar dimensions 
have also presented the values of the DOP close to 1.0. 

Sensitivity of the DOP with the QD stack geom- 
etry parameters: Since the heights of the QDs in the 
experimental TEM images are not very clear (see [lja)), 
so we simulate various geometry configurations of the QD 
stacks and provide the values of the DOP in the table [I] 
This data serve as a measure of the sensitivity of the DOP 
with respect to the QD stack geometry parameters and 
provide a guide to experimentalists to explore the design 
space of such complex multi-million atom systems. From 
the table |TJ as the QD stack height increases in the Li 
— > L3 Lg — > Lg, the value of the degree of polarization 
reduces. The reduction in the value of the DOP is larger 
for the stacks with H=4.0nm as compared to the stacks 
with H=3.5nm. This is due to the fact that the larger 
height of QDs in the stack results in stronger coupling 
between the QDs. This implies a stronger HH-LH inter- 
mixing resulting in larger magnitude of the TMooi-modc. 

The dependence of the DOP on the height (H) of the 
QDs inside the stacks is an unknown factor. The calcu- 
lated values of the DOP in the table Q] shows that the 
DOP becomes very sensitive to the height of the QDs 
inside the stack as the size of the stack grows larger. For 
the systems Li and L3, the increase in the height (H) 
from 3.5nm to 4.0nm results in a small decrease in the 
values of the DOP. However, for a same change in the 
value of H, the DOPno significantly decreases from to 
-0.244 and from 0.445 to -0.45 for the and Lg systems, 
respectively. This implies that an isotropic polarization 
response (DOP ~ 0) can either be achieved from the L 6 
stack with H=3.5nm, or from a stack with fewer number 
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TABLE I: Comparison of experimentally measured and theoretically calculated DOP for various in-plane TE-mode directions 
and QD geometry configurations. Column 1: The multi-layer QD system under study. Column 2: The dimensions of the QDs in 
the stacks: B is base diameter and H is the height of the QD. (V, V) indicates that the QDs are of varying size, base increasing 
by lnm and height increasing by 0.25nm as the size of stack increases in the vertical direction. Column 3-6: The values of the 
DOP calculated from our model. We provide two additional directions for DOP, [100] and [010], for the comparison purpose. 
Column 7-8: The values of DOP computed from experimentally measured TE- and TM-modes presented in the Fig. [ljb-g) or 
taken from T. Inoue et a/P^ 
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of QD layers and H=4.0nm. We therefore propose that 
the polarization response of the QD stacks can be tuned 
by not only increasing the number of QD layer s (a pa - 
rameter tuned in the past experimental studies^ * 5 * 12 * 33 !), 
but also by controlling the height (H) of the individual 
QDs inside the stacks. 

The calculations on the L3 and Lg systems, where the 
size of QDs is increased as the stack size increases, indi- 
cate that such stacks with non-identical QD layers will 
exhibit relatively higher values of the DOP and hence 
will not be suitable for isotropic polarization response. 
The same is true for the case when the height of the QDs 
in the Lg system is increased to 4.5nm such that the ad- 
jacent QD layers touch each other (approaches columnar 
QD limit). Based on the comparison of the calculated 
and the measured values of the DOP for the Lg system 
111 the table |TJ we estimate that the dimensions of the 
QDs inside the stacks are approximately B=20nm and 
H=4nm. 

Conclusions: This article presents a detailed anal- 
ysis of the polarization response of the multi-layer QD 
stacks by PL measurements and through a set of system- 
atic multi-million atom tight binding electronic structure 
calculations. Our theoretical results follow the trends of 
the experimental measurements on quantum dot stacks 
containing single, three, six and nine QD layers and pro- 
vide significant physical insight of the complex physics 
involved by analysing the strain profiles, the band edge 



diagrams, and the wave function plots. The experimen- 
tally measured PL data for the nine quantum dot stack 
reveals a unique property by indicating a significant dif- 
ference in the DOP for the [110] and [110] directions. 
We explain here that this difference is due to the orien- 
tation of the hole wave functions along the [110] direc- 
tion that results in significant reduction of the TE-mode 
along the [110] direction. We suggest that the isotropic 
polarization response from the multi-layer QD stacks is 
due to two factors: (i) the reduction of the TE 110 -mode 
direction and (ii) the increase in the TMnni-mode due 
to enhanced LH-HH intermixing. Our result presented 
in this paper for various geometry configurations serve 
as guidance for the experimentalists to design future QD 
based optical devices. A flip of sign for the DOP in our 
PL measurements and theoretical calculations as the size 
of QD stack increases indicates significant potential to 
achieve polarization insensitive response from multi-layer 
QD stack systems. 
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